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•  Measured  conductivity  of  gadolinia-doped  ceria  by  impedance  spectroscopy. 

•  Described  transport  using  an  equivalent  circuit  including  leads  impedance. 

•  Determined  parameters  from  maxima  and  minima  of  impedance  spectra. 

•  Grain  boundary  arc  can  be  described  accurately  by  a  resistor-capacitor  element. 
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Samples  of  Gd203-doped  Ce02  (GDC)  were  fabricated  by  sintering  of  powder  compacts.  Impedance 
spectra  were  measured  from  400  °C  to  675  °C  in  air  by  electrochemical  impedance  spectroscopy  (EIS). 
Above  ~500  °C,  high  frequency  arc  was  not  semicircular  but  could  be  fitted  with  a  constant  phase 
element  (CPE).  Above  ~625  °C,  high  frequency  arc  could  not  be  resolved  due  to  a  significant  contribution 
from  the  inductive  load.  The  impedance  spectra  were  described  using  a  simple  equivalent  circuit  which 
included  the  leads/instrument  impedance.  The  leads/instrument  impedance  was  measured  over  a  range 
of  frequencies  and  temperatures.  The  high  frequency  part  of  the  impedance  after  subtracting  leads/in¬ 
strument  impedance  could  be  resolved  even  at  the  highest  measurement  temperature  and  was  described 
by  a  semicircle  representative  of  transport  across  grain  boundaries.  From  these  measurements,  grain  and 
grain  boundary  resistivities  were  determined.  The  corresponding  activation  energies  were  0.69  eV  and 
1.11  eV,  respectively.  The  grain  boundary  capacitance  was  nearly  independent  of  temperature.  The 
present  results  show  that  grain  boundary  effects  can  be  described  by  a  resistor  and  a  capacitor.  Relevant 
equivalent  circuit  parameters  were  obtained  from  intercepts,  maxima  and  minima  in  impedance 
diagrams. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electrochemical  Impedance  Spectroscopy  (EIS)  is  a  technique 
which  measures  the  response  of  a  sample  under  an  AC  stimulus  in 
which  the  frequency  is  varied  over  a  wide  range.  EIS  has  wide 
applicability,  and  has  been  used  for  the  study  of  ionic  conductors, 
dielectric  materials,  semiconductors,  solar  cells,  fuel  cells,  batteries 
and  corrosion.  With  improvements  in  equipment  over  the  years, 
four  types  of  techniques  are  typically  used  to  measure  and  interpret 
the  impedance.  They  are:  AC  bridge,  Lissajous  analysis,  phase- 
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sensitive  detection  and  Fourier  analysis.  The  general  approach  in¬ 
volves  the  application  of  a  single  frequency  input  and  the  mea¬ 
surement  of  the  corresponding  current  (/)  flowing  through  and  the 
potential  drop  ( V)  across  the  sample.  From  the  measurement  of  I 
and  V  as  well  as  the  phase  difference  over  a  wide  range  of  fre¬ 
quencies,  the  real  and  the  imaginary  parts  of  the  impedance  as  a 
function  of  frequency  can  be  determined.  In  1969,  Bauerle  in  a 
landmark  paper  in  solid  state  ionics  investigated  transport  through 
yttria-stabilized  zirconia  (YSZ)  samples  with  various  applied  elec¬ 
trodes  and  in  several  atmospheres  using  electrochemical  admit¬ 
tance  spectroscopy  [1],  Bauerle  showed  that  the  measured 
admittance  data  could  be  fitted  to  semicircles  or  depressed  semi¬ 
circles  in  the  admittance  plane.  These  semicircles  at  high  fre¬ 
quencies  correspond  to  transport  properties  through  the  grain 
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interiors  and  across  the  grain  boundaries,  and  at  low  frequencies 
correspond  to  electrode  effects  [1],  Bauerle  used  an  equivalent 
circuit  comprising  a  pure  resistor  and  a  pure  capacitor  in  parallel  to 
describe  the  sample  behavior  at  high  frequencies  which  accurately 
reproduced  the  measured  spectra.  Electrode  effects  were  in  part 
represented  using  Warburg  impedance.  Since  this  original  paper, 
the  use  of  EIS  has  become  commonplace  in  the  study  of  solid  state 
ionics  [2], 

The  resolution  of  grain  and  grain-boundary  effects  usually  re¬ 
quires  a  high  frequency  input,  often  exceeding  100  kHz.  It  is  well 
known  that  at  high  frequencies,  there  can  be  a  significant  contri¬ 
bution  from  an  inductive  load,  the  origin  of  which  lies  in  the  lead 
wires  used  to  connect  to  the  sample  as  well  as  parts  of  the 
instrumentation.  The  effect  of  the  inductive  load  reflects  as  an 
imaginary  part  of  the  impedance  below  the  x-axis  when  the  real 
part  is  plotted  on  the  positive  x-axis  and  negative  of  the  imaginary 
part  is  plotted  on  the  positive  y-axis.  Thus,  in  a  typical  plotting 
scheme,  the  capacitive  effects  appear  above  the  x-axis  and  the 
inductive  effects  appear  below  the  x-axis.  When  both  capacitive 
and  inductive  contributions  are  present  to  a  significant  extent,  their 
effects  are  reflected  in  the  spectra  regardless  of  whether  the  actual 
spectra  lie  above  or  below  the  x-axis.  Many  studies  have  shown  that 
it  is  essential  to  subtract  the  effects  of  the  leads/instrument  before 
data  can  be  adequately  interpreted  [2-4],  In  a  study  of 
(Lai  -xSrxCoi  _yFey03_5)  LSCF  cathodes  using  EIS,  Esquirol  et  al. 
measured  the  leads/instrument  inductance  as  ~  125  nH,  and  the 
corresponding  leads/instrument  impedance  was  subtracted  from 
the  raw  EIS  data  [3],  Samson  et  al.  [4]  also  measured  the  impedance 
of  an  empty  cell  and  subtracted  it  from  the  raw  data  in  their  study 
of  Sr-doped  LaCoC>3  (LSC)  cathodes  [4],  In  many  recent  studies, 
however,  the  approach  used  has  been  to  ignore  that  part  of  the 
impedance  which  lies  below  the  x-axis  and  fit  the  remaining 
portion  above  the  x-axis  (which  is  generally  not  semicircular)  using 
a  CPE  [5—8],  However,  such  a  procedure  does  not  account  for  the 
impedance  of  the  leads/instrument  and  as  such  may  not  represent 
an  accurate  measurement  of  the  sample  impedance.  For  example, 
in  a  study  of  electrolyte  resistivity,  Suzuki  et  al.  did  not  remove  the 
leads/instrument  impedance  [5],  They  noted  that  there  was  a  large 
difference  in  the  measured  electrolyte  resistance  and  that  expected 
from  separate  measurements  on  bulk  samples  [5],  They  attributed 
the  observed  large  difference  in  the  electrolyte  resistance  between 
the  expected  value  and  the  measured  value  to  the  multilayer 
structure  of  their  sample  [5],  However,  the  leads/instrument 
impedance  is  expected  to  contain  both  resistive  and  inductive 
contributions,  which  right-shift  the  arc  leading  to  an  intercept  at  a 
larger  value  (than  the  sample  resistance)  on  the  real  axis  unless  the 
leads/instrument  impedance  has  been  subtracted  out.  Wang  et  al. 
[6]  investigated  Lao.gSro^MnOa-is/Zro^Yo.osCh  porous  electrodes 
by  EIS.  In  their  work  the  high  frequency  intercept  was  equated  to 
the  electrolyte  and  the  lead  resistance.  However,  inductive  effects 
of  the  leads  were  neglected.  Similarly,  Zajac  et  al.  investigated  the 
resistivity  of  doped  ceria  by  EIS  [7,8],  In  their  study,  the  data  below 
the  x-axis  were  cutoff  and  the  remaining  data  were  analyzed  by  a 
CPE  as  the  arcs  clearly  were  not  semicircular  [7,8], 

The  use  of  a  CPE  generally  leads  to  a  better  fit  to  the  data  than 
using  a  capacitor.  This  of  course  is  not  surprising  as  fitting  with  a 
CPE  requires  two  parameters  unlike  a  capacitor,  which  only  needs 
one  parameter.  However,  the  use  of  a  CPE  often  makes  it  difficult  to 
offer  a  unique  physical  interpretation  of  the  data  as  the  fitting  itself 
may  lack  uniqueness.  In  this  work,  we  investigated  the  conductivity 
of  GDC  samples  with  the  objective  being  the  determination  of  grain 
and  grain  boundary  ion  transport  properties.  The  leads/instrument 
impedance  was  separately  measured  and  was  subtracted  from  the 
data  before  fitting.  The  measurements  were  conducted  using  the 
two  terminal  method.  While  many  studies  in  the  past  have 


Fig.  1.  An  SEM  image  of  the  GDC  sample  showing  an  equiaxed  microstructure.  The 
average  grain  size  is  ~  1.1  pm. 

emphasized  the  need  to  subtract  the  leads/instrument  contribution 
[2—4],  the  continual  use  of  procedures  in  which  the  inductive  part 
is  cutoff  and  the  remainder  of  the  data  are  fitted  using  a  CPE  sug¬ 
gests  that  revisiting  this  issue  may  be  of  value.  The  present  study 
also  analyzes  the  data  in  light  of  simple  equivalent  circuits  but  with 
the  inclusion  of  the  leads/instrument  impedance.  The  present  work 
also  attempts  to  extract  relevant  parameters  of  the  equivalent  cir¬ 
cuit  from  intercepts,  maxima  and  minima  in  the  impedance 
diagrams. 

2.  Experimental  procedure 

Pellets  of  Gd2C>3-doped  ceria  (GDC)  powder  of  composition 
80  mol.%  CeC>2— 20  mol.%  GdOi.s  were  die-pressed  and  sintered  in 
air  at  1400  °C  for  2  h.  The  typical  thickness  of  the  pellets  was 
~  1.4  mm.  Pt  paste  was  applied  on  both  sides  of  the  pellets  and  they 
were  fired  at  800  °C  in  air  for  1  h.  Samples  of  two  different  electrode 
geometries  (sizes)  were  used  in  this  study;  one  set  of  samples  with 
electrodes  15.8  mm  in  diameter  and  the  other  with  electrodes 
4.8  mm  in  diameter. 

Impedance  spectra  were  obtained  on  the  samples  in  air  over  a 
range  of  temperatures  between  400  °C  and  675  °C  in  25°  intervals 
using  the  two  terminal  method.  The  impedance  measurements 
were  conducted  using  Solartron  1260  Frequency  Response 
Analyzer  over  a  frequency  (f)  range  between  0.1  Hz  and  1  MHz. 
Impedance  spectra  over  the  same  range  of  temperatures  and  fre¬ 
quencies  were  also  obtained  without  the  sample  (the  empty  cell)  by 
shorting  the  two  lead  wires.  This  corresponds  to  the  inductive- 
resistive  contribution  of  the  lead  wires  and  the  testing  setup. 

3.  Results  and  discussion 

The  density  of  the  sintered  samples  was  measured  to  be 
7.0426  g  ml-1,  which  is  about  97.37%  of  theoretical.  Fig.  1  shows  a 
scanning  electron  micrograph  (SEM)  of  a  polished  and  thermally 
etched  sample  of  a  GDC  sample.  The  microstructure  is  equiaxed  and 
the  average  grain  size  is  ~  1.1  pm.  Fig.  2  shows  an  example  of  how 
the  impedance  was  measured  and  interpreted.  The  corresponding 
equivalent  circuits  are  shown  as  insets.  Also,  an  inset  in  Fig.  2  shows 
schematics  of  the  two  samples  used  with  different  electrode  areas. 
For  simplicity,  the  electrodes  were  also  modeled  using  a  parallel  R— 
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Fig.  2.  Experimental  approach  used  for  impedance  measurements  and  data  interpretation.  Inset  shows  samples  with  two  different  electrode  areas,  (a):  Measured  impedance  of  the 
whole  setup,  including  the  sample  and  the  leads/instrument  at  600  °C;  (b):  Enlarged  view  of  the  high  frequency  region  in  (a);  (c):  Measured  leads/instrument  impedance  by 
shorting  the  leads  without  the  sample;  (d):  Sample  impedance  obtained  by  subtracting  the  leads/instrument  impedance  (c)  from  the  measured  sample-setup  impedance  (a);  (e): 
Enlarged  view  of  the  boxed  part  in  (d).  These  data  are  fitted  to  a  semicircle.  From  the  semicircle,  the  grain  resistance,  the  grain  boundary  resistance  and  the  grain  boundary 
capacitance  are  readily  obtained. 


C  circuit.  However,  it  is  clear  that  the  electrodes  cannot  be  repre¬ 
sented  by  this  simplified  circuit  element.  Other  circuit  elements, 
such  as  Warburg,  and  CPE  will  likely  be  required  to  interpret  the 
low  frequency  data  in  light  of  physically  realistic  mechanisms;  e.g. 
gas  diffusion  through  porous  electrodes,  adsorption,  dissociation, 
charge  transfer.  Thus,  most  of  the  fitting  in  this  manuscript  was 
restricted  to  the  high  frequency  part  of  the  impedance  spectra  even 
though  the  analysis  includes  the  electrodes  modeled  with  a  parallel 
R—C  circuit.  The  low  frequency  part  was  fitted  using  a  CPE  in 
selected  cases  for  completeness  and  to  obtain  a  better  fit  with 
the  data  corresponding  to  the  electrode  effects.  The  measured 


impedance  of  the  entire  setup  including  the  sample  at  600  °C  is 
shown  in  Fig.  2(a).  Fig.  2(b)  shows  an  enlarged  part  of  the  high 
frequency  region  of  Fig.  2(a).  The  impedance  of  the  setup  without 
the  sample  with  the  leads  shorted  is  given  in  Fig.  2(c).  The  leads/ 
instrument  impedance  (Fig.  2(c))  was  then  subtracted  from  the 
measured  total  impedance  (Fig.  2(a)).  The  difference,  which  should 
be  the  sample  impedance  without  the  effects  of  the  leads/instru¬ 
ment  impedance,  is  plotted  in  Fig.  2(d).  Finally,  the  high  frequency 
portion  of  the  impedance  from  Fig.  2(d)  is  plotted  in  Fig.  2(e).  The 
data  points  in  Fig.  2(e)  are  then  fitted  using  a  semicircle.  As  seen  in 
Fig.  2(e),  the  data  can  be  fitted  very  well  with  a  semicircle,  which 
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represents  transport  across  grain  boundaries.  Fitting  with  a  semi¬ 
circle  suggests  that  the  grain  boundary  effects  could  be  adequately 
represented  by  a  resistor  and  a  capacitor  in  parallel.  Also  important 
to  note  is  that  the  data  point  corresponding  to  the  highest  fre¬ 
quency  of  106  Hz  (angular  frequency,  w  =  2it :f =  2n  x  106)  appears 
above  the  x-axis  suggesting  that  the  leads  inductive-resistive  ef¬ 
fects  appear  to  have  been  removed. 

Impedance  data  obtained  on  a  sample  with  electrodes  15.8  mm 
in  diameter  over  a  temperature  range  from  400  °C  to  675  °C  before 
subtracting  the  leads/instrument  impedance  are  plotted  in  Fig.  3(a) 
and  (b).  The  low  frequency  regime,  as  is  well  known,  represents  the 
electrode  effects.  In  this  work  primarily  the  high  frequency  effects 
were  investigated.  At  low  temperatures  (below  525  °C)  the  effects 
of  the  leads/instrument  impedance  are  not  obvious  (although  still 


present)  as  the  entire  measured  spectra  over  the  range  of  fre¬ 
quencies  between  0.1  Hz  and  1  MHz  lie  above  the  x-axis.  The  near 
semicircular  feature  of  the  high  frequency  data  can  be  seen  in 
Fig.  3(a).  For  measurements  at  550  °C  and  above  (Fig.  3(b)),  the  high 
frequency  arc  appears  as  a  depressed  or  a  flattened  semicircle.  At 
the  measurement  temperatures  of  650  °C  and  675  °C,  the  high 
frequency  arc  completely  disappears.  Also,  the  leads/instrument 
inductive  effects  are  clearly  seen  at  higher  temperatures,  as  evi¬ 
denced  by  part  (or  even  most)  of  the  measured  impedance  data 
lying  below  the  x-axis.  The  leads/instrument  impedance 
comprising  inductive  and  resistive  contributions  was  found  to  be 
weakly  dependent  on  temperature,  unlike  the  sample  impedance 
(resistive  part)  which  exhibits  Arrhenius  behavior  and  thus  rapidly 
decreases  with  increasing  temperature.  Thus,  the  effects  of  the 


Fig.  3.  (a)  and  (b)  are  the  as-measured  impedance  spectra  for  the  sample  including  the  leads/instrument  impedance  from  400  °C  to  675  °C;  (c)  and  (d)  are  the  corresponding 
impedance  spectra  after  subtracting  the  leads/instrument  impedance  over  the  same  temperature  range;  (e)  and  (f)  compare  the  plots  with  and  without  the  leads/instrument  effects 
at  650  °C  and  675  °C,  respectively.  Note  that  once  the  leads/instrument  effects  are  subtracted  out,  the  high  frequency  arcs  are  resolved  and  that  they  are  semicircular  in  shape. 
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when  the  time  constants  for  the  three  processes  are  sufficiently  far  apart.  The  impedance  plot  exhibits  3  maxima,  2  minima,  and  4  inflexion  points.  Typically,  the  inflexion  points  are 
located  very  close  to  the  minima.  The  real  part  corresponding  to  w  -*  0  is  given  by  Rj  +  Rg  +  Rgb  +  Re-  The  real  part  corresponding  to  the  high  frequency  intercept  is  equal  to 
Re  z(w*)  =  R|  +  I/RgCg.  (b):  A  schematic  showing  the  expected  impedance  plot  at  an  intermediate  temperature  when  grain,  grain  boundary  and  electrode  effects  can  be  described  by 
parallel  R-C  circuits  and  when  the  time  constants  for  the  three  processes  are  sufficiently  far  apart.  The  impedance  plot  exhibits  2  maxima,  1  minimum,  and  2  inflexion  points.  The 
real  part  corresponding  to  «  ->  0  is  given  by  R[  +  Rg  +  Rgb  +  RE.  The  real  part  corresponding  to  the  high  frequency  intercept  is  given  by  Rez(w')  =  Rj  +  Rg  +  (L  -  R|Cg)/RgbCgb.  (c):  A 
schematic  showing  the  expected  impedance  plot  at  a  high  temperature  when  grain,  grain  boundary  and  electrode  effects  can  be  described  by  parallel  R-C  circuits  and  when  the 
time  constants  for  the  three  processes  are  sufficiently  far  apart.  The  impedance  plot  exhibits  1  maximum.  The  real  part  corresponding  to  <o  -*  0  is  given  by  Ri  +  Rg  +  Rgb  +  RE.  The 
real  part  corresponding  to  the  high  frequency  intercept  is  given  by  Rez(w*)  =  R|  +Rg  +Rgb  +  (L-R|Cg  -  RiCgb)/RECE.  (d):  A  schematic  showing  the  expected  impedance  plot  at 
an  even  higher  temperature  when  grain,  grain  boundary  and  electrode  effects  can  be  described  by  parallel  R-C  circuits  and  when  the  time  constants  for  the  three  processes  are 
sufficiently  far  apart.  The  impedance  plot  exhibits  no  maxima,  no  minima,  and  no  inflexion  points.  The  real  part  corresponding  to  a  -»  0  is  given  by  Ri  +  Rg  +  Rgb  +  RE.  There  is  no  arc 


leads/instrument  impedance  become  more  significant  relative  to 
the  sample  at  higher  temperatures,  which  is  of  course  well  known. 

Fig.  3(c)  and  (d)  shows  the  same  data,  respectively,  from  Fig.  3(a) 
and  (b),  but  from  which  the  leads/instrument  impedance  has  been 
subtracted.  Thus,  the  data  in  Fig.  3(c)  and  (d)  represents  the  sample 
impedance  without  any  leads/instrument  effects.  All  high  frequency 
arcs  now  appear  perfectly  semicircular.  Also,  no  data  points  lie  below 
the  x-axis  showing  that  the  effects  of  the  leads/instrument  have  been 
subtracted  out.  Also  note  that  even  at  the  highest  measurement 
temperature  of  675  °C,  the  high  frequency  semicircular  arc  is  fully 
resolved  and  is  clearly  visible.  Fig.  3(e)  and  (f)  respectively  shows  the 
measured  data  at  650  °C  and  675  °C,  in  which  the  plots  include  data 
with  and  without  leads/instrument  effects.  As  clearly  seen,  the 
semicircle  representing  the  grain  resistance,  the  grain  boundary 
resistance,  and  the  grain  boundary  capacitance  can  be  clearly 
resolved  once  the  leads/instrument  impedance  is  subtracted,  unlike 
the  raw  data  in  which  no  high  frequency  arcs  can  be  seen. 

For  the  high  frequency  data,  such  as  at  1  MHz  noted  in  Fig.  3(e) 
and  (f),  the  change  in  impedance  after  subtracting  the  leads/in¬ 
strument  impedance  consists  of  both  the  real  part  and  the  imagi¬ 
nary  part  corresponding  to  the  leads/instrument  resistance  and  the 
leads/instrument  inductance.  With  decreasing  frequency,  the  effect 
of  the  leads/instrument  inductance  decreases,  but  naturally  the 
resistive  (ohmic)  contribution  remains  the  same.  At  very  low  fre¬ 
quencies,  such  as  at  a  few  Hz  and  below,  the  inductive  part 
essentially  disappears  and  only  the  real  part  from  the  leads/ 


instrument  resistance  remains.  Electrolyte  transport  properties  are 
reflected  typically  in  the  high  frequency  EIS  data,  so  the  leads/in¬ 
strument  impedance  needs  to  be  accounted  for.  Information  con¬ 
cerning  the  electrode  kinetics  is  typically  reflected  in  the  low  to 
intermediate  frequency  EIS  data,  and  thus  one  generally  needs  to 
correct  only  for  the  ohmic  part.  Therefore,  no  matter  whether  the 
electrolyte  or  the  electrodes  are  investigated  in  an  EIS  study,  sub¬ 
traction  of  the  leads/instrument  impedance  is  necessary  when  the 
two  terminal  method  is  used. 

Analysis  of  impedance  spectra  including  electrode,  grain  boundary 
and  grain  effects:  The  most  general  case  but  assuming  simple 
resistor-capacitor  elements  is  examined  in  what  follows.  At  very 
high  frequencies,  the  effects  of  grain  resistance,  Rg,  and  grain 
capacitance,  Cg,  may  become  observable.  This  will  likely  be  the  case 
at  low  temperatures  when  R|Cg  is  large  in  relation  to  (or  compa¬ 
rable  to)  the  leads/instrument  inductance,  L.  The  real  and  the 
imaginary  parts  of  the  impedance  are  given  respectively  by 

Re  2(w)  =  R|  +  7 - — - 7  +  - - — - 7 

('+«!-*?)  (i+^ib) 

+  7 - — - 7  (1) 

(i  +  R!»2ci) 

where  R]  is  the  leads/instrument  resistance  and 
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Imz(w)  M 


JL _ 

[  (i  +Rl<o2q) 

rICe  1 

(l+R2W2C2)J 


(14-^Cg,) 


(2) 


Equation  (2)  shows  that  Im  z(w)  — ►  0  as  w  — ►  0.  The  corre¬ 
sponding  real  part  is  given  by 


time  constants  are  not  sufficiently  different  for  the  three  processes, 
overlapping  spectra  are  expected  and  three  distinct  arcs  may  not  be 
observed.  When  w  =  0,  the  corresponding  real  part  of  the  imped¬ 
ance  is  given  by  Re  z(w  =  0)  =  Ri  +  Rg  +  Rg b  +  Re,  which  is  the  same 
as  Equation  (3).  Assuming  ReCe  »  RgbCgb  »  RgCg,  at  sufficiently 
high  values  of  the  angular  frequency, 


Imz(w)«w 


RgQ  1 

(l  +  R|w2C|)J 


(7) 


Rez(w  —  0)  =  Rj  +  Rg  -t-  Rgb  -t-  Re  (3) 

In  what  follows,  we  will  assume  that  over  some  range  of 
temperatures 


in  which  we  have  substituted  ReCew  >>  RgbCgbw  >>  1.  Thus,  the 
corresponding  w  (which  we  will  denote  w*)  at  which  1m  z(w)  —  0  is 
given  by 


L  <  RgCg  +  R^Cgb  +  ReCe  (4) 

Thus,  there  is  a  range  of  frequencies  over  which  an  arc  exists 
above  the  x-axis  in  the  plot  of  -Im  z(w)  vs.  Re  z(w),  and  at  a  suffi¬ 
ciently  high  angular  frequency,  we  have 


If,  however,  L  >  R2Cg  +  R2bCgb  +  R|Ce,  no  arc  is  observed  above 
the  x-axis.  Such  a  situation  can  occur  in  a  given  system  above  some 
temperature  since  typically,  Rg,  Rgb  and  Re  all  exhibit  an  Arrhenius 
behavior  and  thus  decrease  with  increasing  temperature  while  L  is 
generally  temperature  independent.  Equation  (5)  is  a  cubic  equa¬ 
tion  in  w2  and  should  have  one  real  positive  root,  which  we  will 
denote  by  w*2.  For  angular  frequencies  above  w*  the  Im  z(w)  >  0, 
and  will  appear  below  the  x-axis.  Over  the  angular  frequency  range 
w*  >  w  >  0,  the  Im  z(w)  <  0  and  will  appear  as  an  arc  above  the  x- 
axis.  The  general  shape  of  the  plot  with  —  Im  z(w)  on  the  y-axis  and 
Re  z(w)  on  the  x-axis  is  shown  in  Fig.  4(a). 

The  general  shape  of  the  impedance  spectra  over  the  range 
w*  >  w  >  0  is  characterized  by:  (a)  three  maxima,  (b)  two  minima, 
and  (c)  four  inflexion  points.  These  correspond  to:  (a)  d(-Imz(w))/ 
d(Re  z(w))  =  0,  d2(-Im  z(w))/d(Re  z(w))2  <  0;  (b)  d(-Im  z(w))/ 
d(Re  z(w))  =  0,  d2(-Im  z(w))/d(Re  z(w))2  >  0;  and  (c)  d2(-Im  z(w))/ 
d(Re  z(w))2  —  0.  As  long  as  d(Re  z(w))/dw  is  finite  and  nonzero,  the 
above  conditions  may  also  be  given  as:  (a)  d(-Im  z(w))/dw  =  0, 
d2(-Im  z(w))/dw2  <  0;  (b)  d(-Im  z(w))/dw  =  0,  d2(-Im  z(w))/ 
dw2  >  0;  and  (c)  d2(-Imz(w))/dw2  =  0.  Depending  upon  the  relative 
values  of  the  various  parameters  and  the  range  of  frequencies,  it  is 
possible  to  have:  (i)  Three  maxima,  two  minima,  and  four  inflexion 
points:  (ii)  Two  maxima,  one  minimum,  and  two  inflexion  points; 
(iii)  One  maximum,  and  one  or  two  inflexion  points;  (iv)  One 
maximum;  (v)  No  maxima,  no  minima,  and  no  inflexion  points.  This 
last  situation  corresponds  to  the  case  where  the  entire  spectra 
appear  below  the  x-axis.  In  what  follows,  we  will  first  examine 
some  limiting  cases. 


Case  I.  Low  Temperature:  Suppose  the  temperature  is  sufficiently 
low  such  that 


\l  -  ? - RsCs2  J  =  0  (8) 

L  (l+R2w*2C2)J 

The  w*  is  thus  the  angular  frequency  at  which  capacitive  and 
inductive  parts  of  the  impedance  are  equal  in  magnitude  (but 
opposite  in  sign). 

Thus, 


The  corresponding  real  part  of  the  impedance  is  given  by 

R“(“‘)-  R,+s^  "°> 

At  angular  frequencies  w  >  w*  given  by  Equation  (9),  the 
impedance  spectra  appear  below  the  x-axis.  We  expect  the  general 
expression  for  -Imz(w)  to  be  described  by  Equation  (2)  and  the  plot 
of  — Im  z(w)  vs.  Re  z(w)  given  by  Fig.  4(a).  Such  a  plot  thus  should 
exhibit  three  maxima,  two  minima  and  four  inflexion  points, 
assuming  the  frequency  range  is  wide  enough. 

Case  II.  Intermediate  temperature  range:  At  some  temperature 
(higher  than  for  Case  I),  the  relative  value  of  the  leads/instrument 
inductance  may  be  such  that  R2Cg  <  L  <  R2Cg  +  R2bCgb.  The  point 
w  =  0  gives  the  real  part,  the  same  as  Equation  (3),  and  -lmz(0)  =  0. 
As  the  angular  frequency  is  increased,  there  will  be  a  frequency  at 
which 


“(»•)=“• 


~(1+R>*2CJ,) 


in  which  we  have  substituted  ReCew  »  1  »  RgCgw.  The  corre¬ 
sponding  oj*  is  given  by 


K2hQb  +  K|Cg-L 


(12) 


and  the  corresponding  real  part  of  the  impedance  is  given  by 


L  <  RgCg  (6) 

(and  thus  also  naturally!  <  R2Cg  +  R2bCgb  +  R|Ce).  We  also  assume 
that  ReCe  »  RgbCgb  »  RgCg,  which  are  the  time  constants 
respectively  for  the  electrode,  the  grain  boundary  and  the  grain 
transport  processes.  This  may  lead  to  distinct  arcs  corresponding  to 
electrode,  grain  boundary,  and  grain  effects.  The  corresponding 
observed  behavior  will  be  that  shown  in  Fig.  4(a).  If  however  the 


R“(“") 


R\  +  Rg  + 


(t-RgCg) 

RgbCgb 


(13) 


At  angular  frequencies  higher  than  w*  given  by  Equation  (12), 
the  impedance  spectra  appear  below  the  x-axis.  The  plot 
of  — Im  z(w)  vs.  Re  z(w)  may  likely  exhibit  two  maxima,  one  mini¬ 
mum,  and  two  inflexion  points.  Fig.  4(b)  shows  a  schematic  of  the 
anticipated  E1S  spectra  at  an  intermediate  temperature. 
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We  will  now  examine  the  minimum  in  -Im  z(w)  shown  in 
Fig.  4(b).  The  Im  z(w)  over  the  range  of  frequencies  may  be 
approximately  given  by 


Imz(w)  u,  J\L-  RgCg  -  RgbCgb  -  1  ^  j  (14) 

in  which  we  have  substituted  R|w2C|  «  R2bw2C2b  «  1.  The 
minimum  in  -Im  z(w)  is  given  by 


(21) 


A  plot  of  -Im  z(w)  vs.  Re  z(w)  will  likely  exhibit  one  maximum, 
no  minima,  and  no  inflexion  points.  Fig.  4(c)  shows  a  schematic. 


Case  IV.  An  even  higher  temperature  range:  Suppose  the  leads/in¬ 
strument  inductance,  L,  is  such  that 


l  >  RgCg  +  RgbCgb  +  RlCE 


(22) 


dlmz(te|  _  L-RlCg 
dw  s  8 


which  gives 


a  situation  that  may  occur  at  rather  high  temperatures.  In  such  a 
case,  the  entire  impedance  curve  lies  below  the  x-axis  —  no  arc  is 
seen  above  the  x-axis.  When  o>  =  0,  the  -Im  z(w)  =  0  and  the 
corresponding  real  part  of  the  impedance  is  given  by 


CE  (Rl 


It  is  easily  verified  that  Equation  (16)  corresponds  to  a  minimum 
in  -Imz(w).  Substitution  for  ra2lin  into  the  equation  for  the  real  part 
of  the  impedance  gives  the  real  part  of  the  impedance  corre¬ 
sponding  to  the  minimum  by 


Rez(wmin)  R]  ■  Rg  ■ 


RgbCE(R2Cg  +  R2bCgb-l) 
CE(R2Cg  +  R2bCgb-L)+R2bCg2b 


RECE(R2Cg  +  R2bCgb-L) 
CE(R2Cg+R2bCgb-f.)+R2C2 
Thus  the  arc  length  between  w  =  0  and  o>  —  wmin  is  given  by 


Rez(w  =  0)  -  Rez(w  =  wmin)  = 


CE(R2Cg+R2bCgb-L)+R2bC2b 

_ 

CE(R2Cg+R2bCgb-L)+R2C2 


Note  that  only  if  Rj?Cf » (RgCg+R2bCgb-I)CE  >>R2bC2b,  Equa¬ 
tion  (18)  reduces  to  RE  as  the  arc  length. 


Case  III.  Higher  temperature  range:  At  some  higher  temperature 
(higher  than  for  Case  II),  the  relative  value  of  the  leads/instrument 
inductance  may  be  such  that 


RgCg  +  RgbCgb  <L<  RgCs  +  RgbCgb  +  R|Ce 

The  point  w  =  0  still  gives  the  real  part  the  same  as  Equation  (3) 
and  the  -Im  z(0)  —  0.  As  the  angular  frequency  is  increased,  there 
will  be  a  frequency  at  which 


Imz(w')  -  »[t-KlCs-B|tCit-i+^2j  -  0  (19) 

in  which  we  have  substituted  RgCg w  «  RgbCgbw  <<  1.  The  corre¬ 
sponding  of  is  given  by 


1  \4<=E+R2gbCgb+R2gCg-L 
ReCE  \|  L-R|Cg-R2bCgb 


(20) 


Rez(w  =  0)  =  R{  +  Rg  +  Rgb  +  RE  (23) 

which  is  the  same  as  Equation  (3).  Also,  therefore,  in  this  case 
of  =  0.  Thus,  assuming  data  can  be  realistically  obtained  at  very  low 
frequencies  (w  — »  0),  the  corresponding  intercept  on  the  x-axis 
corresponds  to  the  total  resistance  including  all  components 
regardless  of  the  relative  values  of  the  various  parameters.  Fig.  4(d) 
shows  a  schematic.  Thus,  in  the  event  that  with  increasing  tem¬ 
perature  the  arc  above  the  x-axis  completely  disappears,  the  cor¬ 
responding  intercept  as  «  — >  0  is  still  the  total  resistance. 

Interpretation  of  the  experimental  results  using  the  preceding 
analysis:  In  light  of  the  preceding  analysis,  we  will  now  examine  the 
impedance  spectra  in  Fig.  3(a)  and  (b)  obtained  on  a  GDC  sample 
measured  over  a  temperature  range  from  400  °C  to  675  °C  in  air  and 
over  a  frequency  (f)  range  from  0.1  Hz  to  1  MHz  using  the  two 
terminal  method.  Up  to  525  °C,  almost  the  entire  spectra  are  above 
the  x-axis.  This  means  up  to  525  °C,  the  w*  is  greater  than  the 
maximum  angular  frequency  at  which  measurements  were  made; 
that  is  for  temperatures  up  to  525  °C,  w*  >  2%  x  106.  Above  525  °C, 
part  of  the  measured  spectra  lie  below  the  x-axis,  and  thus  above 
525  °C,  of  <2tz  x  106.  The  spectra  also  show  that  up  to  600  °C,  a 
distinct  high  frequency  arc  can  be  seen.  Thus,  over  the  range  of 
frequencies  investigated,  between  about  425  °C  and  600  °C,  the 
observed  spectra  contain  one  maximum,  one  minimum,  and  two 
inflexion  points.  Also  therefore,  over  the  range  of  experimental 
parameters  (both  the  sample  and  the  testing  parameters),  the  low 
temperature  case  (Case  I,  Fig.  4(a))  is  not  observed.  Over  the  tem¬ 
perature  range  from  400  °C  to  600  °C,  the  observed  spectra 
correspond  to  Case  II  (Fig.  4(b)),  in  which  R2Cg  <  L  <  RgCg  +  RgbCgb. 
Spectra  at  650  °C  and  675  °C  do  not  exhibit  a  high  frequency  arc. 

Over  the  temperature  range  from  400  °C  to  600  °C,  the  high 
frequency  arc  is  clearly  visible.  As  stated  earlier,  the  highest  an¬ 
gular  frequency  at  which  measurements  were  made  was 
Wmax  =  2%  x  106  s-1  in  all  experiments.  Over  the  temperature  range 
from  400  °C  to  450  °C  the  entire  spectra  are  above  the  x-axis.  This 
means  the  corresponding  w*  >  wmax.  The  of  thus  may  be  obtained  by 
extrapolation.  Also,  extrapolation  to  the  x-axis  gives  the  correspond¬ 
ing  real  part  which  is  identified  with  Rez(«*)  —  R|  —  Rg+ 
(L  -  RgCg)/RgbCgb,  namely  Equation  (13).  Over  the  temperature  range 
450  °C-600  °C,  the  high  frequency  arc  can  be  seen  and  at  the  same 
time  some  of  the  data  lie  below  the  x-axis.  In  these  cases,  of  <  wmax. 
The  real  part  of  the  impedance  corresponding  to  of  is  again  given  by 
Equation  (13). 

At  650  °C  and  675  °C,  no  high  frequency  arc  is  observed.  This 
temperature  thus  corresponds  to  the  range  RgCg  +  RgbCgb  <L< 
RgCg  +  RgbCgb  +  R|Ce.  This  corresponds  to  Case  III,  Fig.  4(c).  The  real 
part  of  the  impedance  corresponding  to  w*  is  thus  given  by 
Rez(w’)  -  R,  -  Rg  -  Rgb  +  (L  —  RgCg  -  R2bCgb)/RECE  which  is 


The  corresponding  real  part  is  given  by 
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Table  1 

Fitting  parameters  used  to  fit  simulation  spectra  in  Fig.  3.  Parameters  are  given  only 
up  to  a  temperature  of  600  °C  as  no  minimum  was  observed  at  higher  temperatures. 


T  L  (nH)  Ri  (Q)  Rg(fi)  Cg  Rgb 
fCJ  (F)  fi) 


400  882  3.68  28.4  0  364 

450  874  3.78  10.8  0  107 

500  871  3.88  4.69  0  41 

550  866  4.01  2.37  0  18.4 

600  861  4.1  2.06  0  8.32 


Cgb  tgb  Set  value  xe  (S) 

(nF)  (>iS)  RE(Q)  Ce  (mF) 

34.7  12.6  2.6  k  0.81  2.13 

33.3  3.56  577  1.86  1.07 

30.9  1.26  229  4  0.915 

34.5  0.634  141  5.52  0.776 

38.2  0.318  82.5  7.09  0.585 


T  (°C)  Simulated  Measured  w(kHz)  Simulated  value  ARE (%)  A CE(%) 

Re(Z)  Re  (Z)  RE(fi)  CE 

(A)  OT (mF) 


400  396 

450  122 

500  49.6 

550  24.7 

600  14.5 


413 

129 

51.4 

24.8 


13.4 


0.81  -1.5  -0.19 

1.86  1.22  0.29 

4.01  1.56  0.09 

5.52  -1.2  -0.15 

7.1  1.34  0.04 


Equation  (21 ).  In  view  of  the  fact  that  the  electrode  effects  cannot  be 
adequately  described  by  a  parallel  RE  -  Ce  circuit  element,  the 
applicability  of  Equation  (21 )  to  the  present  data  may  only  be  treated 
as  being  approximate. 

Simulation  and  data  fitting:  Experimental  results  shown  in 
Fig.  3(a)  and  (b)  were  analyzed  using  the  simple  equivalent  circuit 
in  which  electrode,  grain  boundary  and  grain  processes  are 
described  using  parallel  R  -  C  circuits.  The  fitting  of  the  electrode 
effects  by  a  parallel  Re  -  Ce  circuit  is  not  accurate  since  clearly  the 
electrode  arc  is  not  semicircular.  The  primary  focus  of  the  fitting, 
however,  is  the  high  frequency  regime  —  namely  grain  boundary 
and  grain.  Also,  it  was  assumed  that  the  grain  capacitance  Cg  is  very 
small  over  the  range  of  frequencies  and  thus  can  be  set  to  zero 
(actually,  this  means  L  »  R|Cg).  This  assumption,  however,  is  not 
expected  to  be  valid  at  lower  temperatures.  Thus,  the  simulation 
does  not  include  the  low  temperature  Case  I,  which  was  not 
observed  in  the  present  work  over  the  range  of  temperatures  and 
frequencies  investigated. 

Table  1  gives  the  experimental  values  and  the  fitted  values  of  the 
various  parameters  over  a  range  of  temperatures  between  400  °C 
and  700  °C  in  50°  intervals.  Experimental  data  in  Fig.  3  correspond 
to  one  minimum  below  600  °C  (the  schematic  corresponding  to 
Fig.  4(b)).  This  minimum  corresponds  to  Case  II,  the  intermediate 
temperature  case.  The  corresponding  wjCn  is  given  by  Equation 
(16).  Table  1  shows  the  values  of  L,  R|,  Rg,  Rgb,  Cgb,  Re  and  Ce  (with  Cg 
assumed  to  be  zero).  Note  that  the  L  is  essentially  independent  of 
temperature  (it  very  slightly  decreases  with  increasing  tempera¬ 
ture)  and  R]  slightly  increases  with  increasing  temperature.  Table  1 
also  lists  wmin  over  the  temperature  range  from  400  °C  to  600  °C.  At 
650  °C  and  above,  which  corresponds  to  Case  III,  no  high  frequency 
arc  is  observed,  and  wjCn  <  0,  which  means  no  minimum  exists 
(that  is  at  and  above  650  °C,  L  >  RgbCgb).  The  values  of  the  various 
parameters  in  Table  1  indeed  show  that  at  and  above  650  °C, 
L  >  Rgb Cgb-  The  estimated  values  of  Re  and  Ce  from  the  simulation 
are  very  close  to  the  set  values  of  Re  and  Ce. 

Fig.  5(a)  through  (f)  is  the  simulated  EIS  spectra  over  the  tem¬ 
perature  range  400  °C-650  °C  in  50°  intervals  and  over  the  fre¬ 
quency  range  from  0.1  Hz  to  106  Hz.  Also  shown  in  each  of  the 
figures  (from  Fig.  5(a)  through  (e))  is  an  inset  corresponding  to  the 
minimum  in  -Im  z(w)  (corresponding  to  wmjn).  From  each  of  the 
minima,  the  corresponding  Re  z(«mjn)  was  estimated,  also  shown  in 
the  insets.  The  experimental  values  of  Re  z(wmin)  were  determined 
from  the  measured  spectra  given  in  Fig.  3(a)  and  (b).  These  values 
are  also  listed  in  Table  1.  The  agreement  between  the  two  is  very 
good.  In  Fig.  5(f),  which  shows  spectra  at  650  °C,  no  high  frequency 


arc  is  observed.  The  inset  shows  enlarged  view  of  the  data  near  the 
intersection  with  the  x-axis. 

Table  1  also  lists  the  values  of  Rg  and  Rgb  as  a  function  of  tem¬ 
perature.  The  Rg  varies  between  ~28.8  Q  at  400  °C  and  ~2.04  Q  at 
575  °C.  Over  the  same  range  of  temperatures,  the  Rgb  varies  be¬ 
tween  ~364  Q  and  11.94  fi.  These  results  show  that  even  at  575  °C, 
the  Rgb  is  greater  than  Rg.  Above  about  600  °C,  the  Rg  appears  to  be 
independent  of  temperature.  But  this  observation  is  related  to 
inaccuracies  in  the  measurements  (errors  involved  in  the  removal 
of  the  leads  impedance).  It  is  the  expectation  that  both  Rg  and  Rgb 
will  continue  to  decrease  with  increasing  temperature.  In  order  to 
improve  the  accuracy  of  measurements,  one  would  need  to  use 
samples  of  larger  total  resistance  (smaller  electrode  area  and/or 
thicker  samples)  compared  to  the  ones  used  in  the  present  work. 
The  Rgb  continues  to  decrease  with  increasing  temperature.  There 
likely  are  greater  uncertainties  in  both  values  above  600  °C.  The 
observation,  however,  that  even  at  a  temperature  as  high  as  650  °C 
the  Rgb  is  comparable  to  (or  greater  than)  Rg  implies  that  the 
absence  of  a  high  frequency  arc  in  the  measured  spectra  (when  not 
corrected  for  the  leads/instrument  impedance)  cannot  be  inter¬ 
preted  to  mean  that  the  grain  boundary  resistance  is  negligible 
compared  to  the  grain  resistance. 

Table  1  shows  that  the  Cgb  is  essentially  independent  of  tem¬ 
perature  over  the  range  400  °C-600  °C,  and  seems  to  increase 
slightly  above  600  °C  This  apparent  increase  in  Cgb  at  higher 
temperatures  is  again  related  to  the  inaccuracies  in  the  measure¬ 
ments  and  samples  of  appropriate  geometries  will  be  required  to 
obtain  accurate  values  of  all  three  parameters,  namely  Rg,  Rgb  and 
Cgb-  The  Cgb  is  expected  to  be  essentially  independent  of  tempera¬ 
ture  over  a  very  wide  range  since  impurity  and  dopant  profiles  near 
the  grain  boundaries  and  the  associated  space  charge  effects  are 
practically  frozen  in  from  the  sintering  temperature  [9-12], 

Electrode  parameters :  As  stated  earlier,  the  electrode  effects 
cannot  be  accurately  described  by  a  parallel  R—C  circuit  element 
given  the  nature  of  many  complex  processes  that  occur  at  the 
electrodes,  such  as  gas  phase  diffusion  through  porous  electrodes, 
adsorption  and  dissociation  of  02  molecules,  electron  transfer 
corresponding  to  the  oxygen  reduction  reaction  (ORR)  as  well  as 
the  oxygen  ion  oxidation  reaction  (the  reverse  of  ORR),  and  the 
incorporation/extraction  of  O2-  at  the  electrode/electrolyte  inter¬ 
face  (typically  at  three  phase  boundaries,  TPBs).  However,  if  elec¬ 
trode  effects  can  be  adequately  described  in  any  given  system  by  a 
parallel  R—C  circuit,  the  present  work  shows  that  it  may  be  possible 
to  obtain  both  Re  and  Ce  from  the  minimum  in  —  Im  z(«)  when  it  is 
impractical  to  make  measurements  at  ultra-low  angular  fre¬ 
quencies.  Table  1  compares  the  estimated  values  of  Re  and  CE  from 
the  minima  in  -Im  z(w)  with  those  used  in  the  simulations.  It  is 
observed  that  the  error  in  Re  is  typically  less  than  2%  and  that  in  Ce 
is  much  lower.  In  the  present  work,  measurements  were  made  at 
frequencies  >0.1  Hz,  which  precluded  the  estimation  of  the  Re  from 
the  arc  length.  Thus,  the  values  of  electrode  parameters,  namely  RE 
and  Ce,  were  obtained  from  the  minimum,  especially  at  lower 
temperatures.  An  examination  of  Fig.  3(a)  and  (b)  shows  that  if  the 
low  frequency  arcs  are  extrapolated  to  zero  angular  frequency  (to 
the  x-axis),  the  values  obtained  for  Re  are  on  the  same  order  of 
magnitude  as  selected  for  the  simulation  and  the  ones  estimated 
from  the  minima  in  -Im  z(w).  Although  these  estimates  are  not 
expected  to  be  accurate,  some  conclusions  may  nevertheless  be 
drawn  from  the  values  obtained  for  Re  and  Ce. 

Over  the  temperature  range  from  400  °C  to  600  °C,  the  Re  varies 
between  2620  Q  and  82.45  Q,  or  by  a  factor  of  ~  31,  clearly  showing 
that  the  overall  electrode  reaction  exhibits  a  strong  temperature 
dependence.  It  is  well  known  that  the  reaction  of  charge  transfer, 
described  by  such  models  as  the  Butler— Volmer  equation,  is  usually 
thermally  activated,  and  a  strong  temperature  dependence  is  to  be 
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expected,  as  observed  here.  The  other  temperature-dependent 
processes  which  affect  Re  include  gas  phase  adsorption  and  disso¬ 
ciation.  Indeed,  in  the  original  landmark  paper  by  Bauerle  [1],  the 
resistance  describing  the  electrode  reaction  exhibited  a  thermally 
activated  behavior.  Over  the  same  range  of  temperatures,  the 
electrode  capacitance,  Ce  varies  from  8.11  x  10-4  f  at  400  °C  to 
7.09  x  10-3  f  at  600  °C.  That  is  the  Ce  varies  by  ~9  times.  Such  a 
large  variation  in  Ce  is  consistent  with  the  nature  of  processes  that 
occur  at  the  electrodes.  One  of  the  important  processes  that  occurs 
at  the  electrolyte/electrode  interface  is  that  of  O2  adsorption  and  its 
dissociation.  This  determines  the  surface  coverage  of  adsorbed  02 
(or  O)  and  thus  would  dictate  the  overall  charge  transfer  reaction. 
The  temperature  dependence  of  the  adsorption  and  dissociation 
processes  would  thus  be  expected  to  affect  the  Ce  resulting  in  a 
temperature  dependence  of  Ce.  Thus,  unlike  the  grain  boundary 


capacitance,  Cgb,  which  is  expected  to  be  essentially  temperature- 
independent,  the  electrode  capacitance,  Ce,  would  likely  exhibit 
temperature  dependence.  As  stated  earlier,  the  description  of 
electrode  effects  by  a  parallel  resistor-capacitor  circuit  element  is 
an  over  simplification.  In  principle,  however,  it  should  be  possible 
to  develop  similar  analysis  as  given  here  but  with  a  more  repre¬ 
sentative  circuit  element  for  the  electrodes  (e.g.  CPE).  If  this  can  be 
done,  it  may  be  possible  to  estimate  the  electrode  parameters  from 
the  analysis  of  the  minima  in  -Im  Z(oj)  vs.  to;  that  is  from  wmin. 

The  analysis  given  here  assumes  that  the  time  constants  for  the 
three  processes,  namely,  grain,  grain  boundary  and  electrode,  are 
sufficiently  far  apart  so  that  impedance  spectra  should  exhibit 
distinct  arcs.  In  the  simulation,  the  grain  capacitance,  Cg  was  set  to 
zero;  that  is,  its  time  constant  was  set  to  zero.  From  the  values  of  Re 
and  Ce,  the  estimated  time  constant  for  the  electrode  process, 
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Fig.  6.  Comparison  between  leads/instrument  impedance  subtraction  and  inductive  cut-off,  as  well  as  between  capacitance-fitting  and  CPE-fitting,  (a):  The  high  frequency  regime 
of  the  measured  impedance  with  large  Pt  electrodes  at  600  °C.  (b):  The  high  frequency  arc  from  (a)  after  cutting  off  the  data  below  the  x-axis.  The  arc  is  not  semicircular  and  was 
fitted  with  a  CPE.  (c):  The  high  frequency  arc  from  (a)  after  subtracting  the  leads/instrument  impedance,  (d):  The  high  frequency  arc  from  (c).  The  arc  could  be  accurately  fitted  as  a 


te  =  ReCe,  varies  between  ~2.12  s  at  400  °C  and  —0.59  s  at  600  °C 
Similarly,  from  the  values  of  Rgb  and  Cgb,  the  estimated  time  con¬ 
stant  for  the  grain  boundary  process,  igb  =  RgbCgb,  varies  between 
~  1.26  x  10"5  s  at  400  °C  and  -3.18  x  10"7  s  at  600  °C.  Clearly,  the 
time  constants  for  the  electrode  and  the  grain  boundary  processes 
are  sufficiently  far  apart  which  should  lead  to  distinct  arcs.  This  also 
means  that  the  disappearance  of  the  high  frequency  arc  at  high 
temperatures  is  not  related  to  any  overlapping  spectra;  but  is  rather 
related  to  a  large  contribution  from  the  leads/instrument  inductive 
effects. 

The  typical  electrode  thickness  was  about  25  microns.  Also,  the 
electrodes  exhibited  a  sufficiently  high  level  of  porosity  (typically 
-30%  based  on  prior  work).  The  approximate  estimate  of  the  time 
constant  (not  strictly  a  time  constant  as  the  corresponding  differ¬ 
ential  equation  does  not  lead  to  a  single  time  constant)  for  gas 
diffusion  through  the  porous  electrodes  is  estimated  to  be  less  than 
0.01  s.  The  origin  of  the  estimated  long  time  constants  for  the 
electrode  processes  probably  lies  in  sluggish  electrode  reaction  of 
charge  transfer  (including  adsorption,  dissociation  and  charge 
transfer). 

A  comparison  between  leads  I  instrument  impedance  subtraction 
and  inductive  cutoff:  The  measured  impedance  spectra  at  600  °C 
given  in  Fig.  2(b)  are  reproduced  in  Fig.  6(a)  along  with  the 
equivalent  circuit  used  for  fitting  (the  high  frequency  regime).  The 
inductive  effects  at  high  frequencies  are  clearly  seen  in  the  figure.  In 
Fig.  6(b),  the  data  below  the  x-axis  are  cutoff  and  an  arc  is  fit  to  the 


data  using  a  CPE.  Note  that  the  arc  is  not  semicircular,  as  to  be 
expected  due  to  the  presence  of  inductive  effects.  In  accord  with  the 
equivalent  circuit  analysis  given,  the  high  frequency  intercept  is 
given  approximately  by  Ri  +  Rg  +  L/RgbCgb.  The  magnitude  of  the 
intercept  from  Fig.  6(b)  is  —  8.45  O.  The  magnitude  of  the  intercept 
estimated  using  values  from  Table  1  is  —8.87  Q  showing  good 
agreement.  The  low  frequency  intercept  from  Fig.  6(b)  is  - 13.95  Q. 
The  value  of  the  low  frequency  intercept  estimated  from  Table  1  is 
- 14.48  Q,  again  showing  good  agreement.  The  arc  length  from  the 
extrapolation  of  the  data  is  —  5.5  Q,  while  that  based  on  values  from 
Table  1  is  —5.61  Q,  which  shows  very  good  agreement. 

Fig.  6(c)  shows  the  high  frequency  data  from  which  the 
measured  leads/instrument  impedance  was  subtracted.  As  seen  in 
the  figure,  the  entire  spectra  now  lie  above  the  x-axis.  In  Fig.  6(d) 
the  high  frequency  part  of  the  arc  from  Fig.  6(c)  is  reproduced,  and 
to  which  a  semicircular  arc  is  fitted.  Note  that  the  arc  can  be  fitted  to 
a  semicircle  quite  well.  The  high  frequency  intercept  in  Fig.  6(d)  is 
identified  with  Rg.  From  Fig.  6(d),  the  estimated  value  of  Rg  is 
—  2.15  Q.  From  Table  1,  the  value  of  Rg  is  —2.06  Q,  showing  very 
good  agreement.  The  low  frequency  intercept  in  Fig.  6(d)  is  iden¬ 
tified  with  Rg  +  Rgb,  whose  magnitude  is  — 10.3  Q.  From  Table  1,  the 
estimated  value  of  Rg  +  Rgb  is  - 10.38  Q,  again  showing  very  good 
agreement.  The  preceding  shows  that  after  subtraction  of  the  leads/ 
instrument  impedance,  the  high  frequency  data,  which  represent 
grain  and  grain  boundary  effects,  can  be  well  represented  by  a 
parallel  Rgb  -  Cgb  circuit  element. 
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CPE  fits  to  low  frequency  data :  Fig.  3(e)  (650  °C)  and  (f)  (675  °C) 
are  reproduced  respectively  as  Fig.  7(a)  and  (b).  The  points  in  black 
correspond  to  the  as-measured  data.  The  points  in  red  correspond 
to  the  data  after  subtracting  the  leads/instrument  impedance.  The 
high  frequency  semicircular  arcs  are  well  resolved  once  the  leads/ 
instrument  impedance  is  removed.  Fig.  7(c)  and  (d)  shows  the 
corresponding  simulated  spectra  in  which  the  low  frequency 
electrode  effects  are  described  by  a  CPE.  The  corresponding  pa¬ 
rameters  are  given  in  Table  2.  As  seen,  the  simulated  spectra  are  in 
good  agreement  with  the  observed  ones.  However,  as  stated  earlier, 
two  parameters  are  needed  for  a  CPE  fit  unlike  a  capacitor  which 
requires  only  one.  While  the  fits  with  a  CPE  can  be  well  matched 
with  the  experimental  data,  physical  interpretation  is  often  elusive. 
It  would  appear  that  equivalent  circuit  elements  that  represent 
relevant  physical/chemical  phenomena  are  necessary.  Results  of 
fitting  with  a  CPE  are  given  only  to  demonstrate  that  electrode 
effects  can  be  adequately  fitted  with  a  parallel  CPE-resistor  circuit. 

Estimation  of  grain  and  grain  boundary  parameters  using  the 
measured  spectra:  In  order  to  determine  Rg,  Rgb,  and  Cgb,  the  leads/ 
instrument  impedance  was  subtracted  from  the  measured  imped¬ 
ance  spectra.  Each  of  the  high  frequency  arcs,  once  the  leads/in¬ 
strument  impedance  was  subtracted  out  (Fig.  3),  could  be 
adequately  described  by  a  semicircle.  In  all  cases,  the  semicircular 
arcs  were  extrapolated  to  the  x-axis.  The  high  frequency  intercept 
of  the  semicircle  corresponds  to  the  grain  resistance,  Rg,  and  the 
low  frequency  intercept  corresponds  to  the  grain  +  grain  boundary 
resistance,  Rg  +  %,  [  1  ].  The  corresponding  Arrhenius  plots,  namely, 
ln(Rg/T)  vs.  1000/T  and  ln(Rgb/T)  vs.  1000/T,  are  shown  in  Fig.  8.  In 
Fig.  8(a)  and  (b)  first  the  raw  data  are  shown  on  Arrhenius  plots, 


Table  2 

Parameters  used  for  fitting  data  at  650  °C  and  675  °C  from  Fig.  3  using  a  CPE  for  the 
low  frequency  arc.  The  corresponding  CPE  fits  to  the  data  are  shown  in  Fig.  7(a)  and 
W. 

T  L  Ri  Rg  Rgb  Re  Cgb  CPE-T  CPE-P 

650  8.58E-7  4.2  2.16  3.4  35.99  6.62E-8  9.4E-3  0.48 

675  8.57E-7  4.25  2.28  1.88  22.24  1.2E-7  0.01  0.51 


from  which  the  leads/instrument  impedance  was  not  subtracted 
out.  The  data  exhibit  a  linear  behavior  on  the  Arrhenius  plots  only 
up  to  a  temperature  of  ~500  °C.  Above  500  °C,  the  high  frequency 
arc  is  not  clearly  resolved  in  the  impedance  spectra.  This  reflects  as 
an  apparent  deviation  from  the  linear  behavior  on  the  Arrhenius 
plots  and  the  corresponding  error  is  ~  8  Q.  It  is  to  be  noted  that  the 
non-removal  of  the  leads/instrument  impedance  means  that 
despite  Arrhenius-like  behavior  below  500  °C,  the  measured  acti¬ 
vation  energies  likely  contain  errors  related  to  the  non-removal  of 
the  leads/instrument  impedance. 

Fig.  8(c)  and  (d)  shows  the  same  data  after  subtracting  the  leads/ 
instrument  impedance.  The  linear  behavior  on  Arrhenius  plots  is 
now  preserved  up  to  a  temperature  of  600  °C  and  the  corre¬ 
sponding  error  level  is  reduced  to  ~  2  Q.  The  data  shown  in  Fig.  8(a) 
through  (d)  were  obtained  on  a  sample  with  large  platinum  elec¬ 
trodes  (15.8  mm  diameter)  and  with  platinum  leads.  Finally, 
Fig.  8(e)  and  (f)  shows  data  obtained  using  small  platinum  elec¬ 
trodes  (4.8  mm  diameter)  and  with  silver  leads,  after  subtracting 
out  the  leads/instrument  impedance.  The  linear  behavior  is  now 
observed  up  to  ~675  °C  and  the  error  level  is  further  reduced  to 
~6  Q.  It  is  to  be  noted  that  on  a  relative  basis  an  error  of  6  Q  with 


Fig.  7.  EIS  spectra  at  650  °C  and  675  °C.  (a)  650  °C  measured  spectra  and  after  subtracting  the  leads/instrument  impedance,  (b)  675  °C  measured  spectra  and  after  subtracting  the 
leads/instrument  impedance,  (c)  650  °C  simulated  spectra  (with  a  CPE)  with  and  without  the  leads/instrument  impedance,  (d)  650  °C  simulated  spectra  (with  a  CPE)  with  and 
without  the  leads/instrument  impedance. 
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Fig.  8.  Arrhenius  plots  corresponding  to  ln(Rg/T)  vs.  1000/r  (a,  c,  e)  and  ln(Rgb/T)  vs.  1000/T  (b,  d,  f)  where  Rg  is  the  grain  resistance  and  Rgb  is  the  grain  boundary  resistance,  (a)  and 
(b)  correspond  to  the  intercepts  obtained  with  large  Pt  electrodes  and  Pt  leads  without  subtracting  the  leads/instrument  impedance,  (c)  and  (d)  correspond  to  the  intercepts 
obtained  using  data  from  (a)  and  (b)  but  after  subtracting  the  leads/instrument  impedance,  (e)  and  (f)  correspond  to  the  intercepts  obtained  with  small  Pt  electrodes  and  Ag  leads 
after  subtracting  the  leads/instrument  impedance. 


smaller  electrodes  is  actually  lower  than  an  error  of  2  0  with  larger 
electrodes.  The  present  results  thus  show  that  by  a  suitable  choice 
of  a  sample  geometry  and  subtracting  the  leads/instrument 
impedance,  the  high  frequency  arc  could  be  recovered  and  could  be 
fitted  well  to  a  semicircle  consistent  with  the  description  of  the 
grain  boundary  effects  using  a  resistor  and  a  capacitor,  as  originally 
shown  by  Bauerle  in  his  landmark  paper  [1],  It  should  be  further 
possible  to  choose  a  sample  of  an  appropriate  geometry  to  increase 
its  resistance  relative  to  the  leads/instrument  impedance  and 
extend  the  linear  regime  of  the  Arrhenius  plots  to  even  higher 
temperatures. 

The  data  from  Fig.  8(a),  (c)  and  (e)  can  be  fitted  to  [1] 


-  «STexp(^)  (24) 

for  describing  the  temperature  dependence  of  the  grain  resistance 
with  R°  as  the  pre-exponent  (units  of  Q  K-1)  and  Qg  as  the  activa¬ 
tion  energy.  Similarly,  data  from  Fig.  8(b),  (d)  and  (f)  can  be  fitted  to 


for  describing  the  temperature  dependence  of  the  grain  boundary 
resistance  with  R“b  as  the  pre-exponent  (units  of  0  K-1)  and  Qgb  as 
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Table  3 

Activation  energies  for  transport  across  grain  and  grain  boundaries  of  the  GDC 
samples  used  in  the  present  study.  Measured  activation  energies  are  compared  with 
literature  data  (Ref.  [13]). 


Raw 

L-Pt 

L-Ag 

S-Ag 

Literature  [13] 


Qg(eV) 


Qgb(eV) 


0.67 

0.73 


0.78 


the  activation  energy.  The  corresponding  activation  energies  are 
listed  in  Table  3.  The  data  obtained  with  small  electrodes  and  silver 
wires  correspond  to  the  largest  linear  range  on  the  Arrhenius  plots. 
For  the  sample  with  large  electrodes  (Fig.  8(c)  and  (d)),  the  corre¬ 
sponding  values  were  estimated  to  be:  R°  =  1.2275  x  10-7  Q  K-1  and 
Qg  —  0.73  eV  for  the  grain  resistance  and  R°b  =  1.3684  x  10-8  Q  K-1 
and  Qgb  =  1.01  eV  for  the  grain  boundary  resistance.  For  the  sample 
with  small  electrodes  (Fig.  8(e)  and  (f)),  the  corresponding  values 
were  estimated  to  be:  R°  =  2.0148  x  10  6  Q  K-1  and  Qg  =  0.69  eV  for 
the  grain  resistance  and  R°b  =  1.5741  x  10-8  Q  K-1  and  Qgb  =  1.11  eV 
for  the  grain  boundary  resistance.  The  measured  activation  energies 
are  compared  with  literature  values  in  Table  3  [13],  The  values  of  the 
activation  energies  determined  in  the  present  work  from  samples 
with  small  electrodes,  namely  Qg  =  0.69  eV  and  Qgb  =  1.11  eV,  are 
deemed  more  accurate  than  those  measured  using  larger  electrodes. 

From  Equations  (24)  and  (25),  one  may  estimate  grain  and  grain 
boundary  resistances  of  the  samples  at  other  temperatures.  At 
800  °C,  the  estimated  values  of  the  grain  resistance  (small  elec¬ 
trodes)  is  3.764  Q  and  that  of  the  grain  boundary  resistance  is 
2.762  Q.  Thus,  approximately  ~58%  of  the  resistance  at  800  °C  is 
attributed  to  the  grain  resistance  and  ~42%  of  the  resistance  is 
attributed  to  the  grain  boundary  resistance.  At  lower  temperatures, 
the  contribution  of  the  grain  boundary  resistance  will  be  even 
greater.  Thus,  the  present  work  shows  that  in  order  to  lower  the 
electrolyte  resistance,  it  is  necessary  to  reduce  the  grain  boundary 
contribution.  The  main  contributors  to  grain  boundary  resistance 
are  the  grain  size  (the  smaller  the  grain  size,  the  higher  is  the  grain 
boundary  contribution  to  the  total  resistance),  possible  impurities 
at  grain  boundaries  and  space  charge  effects.  Note  that  the  sample 
grain  size  in  the  present  work  was  rather  small,  ~1.1  pm. 

From  the  high  frequency  semicircle,  the  grain  boundary  capac¬ 
itance  can  be  obtained  from  the  maximum  in  -Imz(w).  This  is  given 
by  [l] 


WmRgb 


(26) 


where  wm  is  the  angular  frequency  at  which  -1m  z(w)  is  a 
maximum. 

The  Cgb  determined  using  the  data  obtained  with  Pt  leads  with 
large  Pt  electrodes  and  Ag  leads  with  small  Pt  electrodes  are  plotted 
in  Fig.  9  vs.  the  measurement  temperature.  The  data  obtained  with 
large  electrodes  show  an  increase  with  increasing  temperature  above 
about  600  °C.  This,  however,  is  attributed  to  errors  associated  with 
the  incomplete  removal  of  the  leads/instrument  impedance  which 
results  in  significant  errors  in  resolving  the  spectra  at  higher  tem¬ 
peratures  for  the  aforementioned  reasons.  The  data  obtained  with 
small  electrodes,  however,  show  that  the  measured  grain  boundary 
capacitance  is  essentially  independent  of  temperature  over  the  range 
of  temperatures  over  which  data  were  obtained.  This  is  consistent 
with  expectations  in  that  the  grain  boundary  capacitance  is  likely 
determined  by  the  grain  boundary  structure  and  chemistry  [9—12], 
which  is  virtually  frozen  in  at  low  temperatures  over  which  the 
measurements  were  made.  Figs.  8  and  9  respectively  give  absolute 


T  (°C) 


Fig.  9.  Grain  boundary  capacitance,  Cgb  obtained  by  fitting  to  the  high  frequency  semi¬ 
circle.  L_Pt  corresponds  to  the  sample  with  large  electrodes  and  measured  using  Pt  leads; 
S_Ag  corresponds  to  the  sample  with  small  Pt  electrodes  and  measured  using  Ag  leads. 


values  of  Rg  &  Rg b  and  Cgb  instead  of  area  specific  values.  The  differ¬ 
ence  in  the  capacitance  values  for  electrodes  of  smaller  and  larger 
diameters  are  in  part  due  to  the  sizes  of  the  electrodes.  There  is  also 
an  edge  effect,  because  of  which  neither  the  capacitance  nor  the 
resistance  values  scale  exactly  with  the  form  factors  (ratio  of  the 
electrode  area  to  the  sample  thickness)  of  the  electrodes. 

4.  Summary 

EIS  measurements  were  made  on  GDC  disc-shaped  samples 
with  Pt  electrodes  of  two  different  sizes;  4.8  mm  in  diameter  and 
15.8  mm  diameter.  Measurements  were  made  over  a  range  of 
temperatures  between  400  °C  and  675  °C  in  25°  intervals.  Above 
~500  °C,  the  high  frequency  spectra  of  the  samples  appeared  to  be 
flattened  semicircles.  At  650  °C  and  675  °C,  no  high  frequency  arcs 
were  observed.  Also,  a  significant  portion  of  the  measured 
impedance  at  high  frequencies  was  below  the  x-axis  consistent 
with  a  significant  contribution  from  the  leads/instrument  induc¬ 
tance.  The  leads/instrument  impedance  was  separately  measured 
by  shorting  the  lead  wires  without  the  sample.  Once  the  lead 
impedance  was  subtracted  from  the  measured  total 
(sample  +  leads/instrument)  impedance,  the  high  frequency  arcs 
could  be  recovered  at  the  highest  measurement  temperature  of 
675  °C.  Also,  the  arcs  could  be  fitted  accurately  with  a  semicircle  in 
each  case,  consistent  with  the  corresponding  equivalent  circuit 
being  a  resistor  and  a  capacitor  in  parallel  corresponding  to  trans¬ 
port  across  the  grain  boundaries.  From  the  measurements,  after 
subtracting  the  leads/instrument  impedance,  the  grain  resistance, 
Rg,  the  grain  boundary  resistance,  Rgb,  and  the  grain  boundary 
capacitance,  Cgb,  could  be  resolved.  The  grain  resistance  and  the 
grain  boundary  resistance  exhibited  thermally  activated  behavior 
with  the  corresponding  activation  energies  of  Qg  =  0.69  eV  and 
Qgb  =  1.11  eV,  respectively.  The  grain  boundary  capacitance  was 
essentially  independent  of  temperature.  The  present  work  un¬ 
derscores  the  importance  of  subtracting  the  leads/instrument 
impedance,  in  accord  with  some  earlier  studies  [2—4],  Once  the 
leads/instrument  impedance  is  subtracted  out,  the  high  frequency 
arc  could  be  accurately  fitted  with  a  capacitor  (one  parameter)  and 
a  resistor,  and  it  was  not  necessary  to  use  a  CPE  (which  requires  two 
fitting  parameters). 

The  as-measured  spectra  were  also  fitted  to  a  simple  equivalent 
circuit  in  which  grain,  grain  boundary,  and  electrode  effects  were 
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represented  by  parallel  R—C  circuit  elements.  The  representation  of 
the  electrode  effects  by  a  parallel  R—C  circuit  is  an  over  simplification. 
Nevertheless,  many  salient  features  of  the  observed  spectra  could  be 
deduced  by  analyzing  maxima  of  the  -Im  z(«)  vs.  Re  z(w)  plots  as 
well  as  minima  of  the  -1m  z(w)  vs.  Re  z(w)  plots.  It  was  shown  that 
from  the  minima  in  -Im  z(w)  vs.  Re  z(«)  plots,  one  can  obtain  pa¬ 
rameters  representing  the  equivalent  circuit.  Thus,  in  cases  in  which 
the  experimental  range  is  not  sufficient  to  estimate  the  parameters 
from  the  maxima  and  the  intercepts,  such  as  for  example  Re  and  Ce, 
approximate  estimates  may  be  obtained  from  the  minima.  The  pre¬ 
sent  work  also  shows  that  an  apparent  disappearance  of  the  grain 
boundary  arc  in  measured  EIS  spectra  at  high  frequencies  is  likely  the 
result  of  not  having  removed  the  leads/instrument  impedance  and 
cannot  always  be  interpreted  to  mean  that  the  grain  boundary 
resistance  is  negligible  in  relation  to  the  grain  resistance. 
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